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In radiation penetration theory, infinite medium flux distributions have for some years been 
calculated using biorthogonal functions called U k n {z). In this paper the spaces spanned by these func- 
tions, transformation kernels, generating functions, recursion relations, asymptotic trends for large 
w, many expansions, and relations with well-known orthogonal polynomials are worked out. 
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1. Introduction 

The biorthogonal (BON) system discussed in this paper was encountered in 1950 during at- 
tempts to construct spatial distributions of the gamma-ray flux from the values of the first few 
even moments. The desired distributions appeared to depend on |z|, rather than z 2 , so that it did 
not seem reasonable to represent the distributions with polynomials which utilized powers of z 2 , 
as would have been the case for an orthogonal polynomial representation. In discussing this basic 
asymmetry with Spencer, U. Fano suggested that perhaps one should approximate with functions 
whose Laplace transforms were even functions of the transform parameter, but otherwise with 
the right kind of singularity, properly located. 

This suggestion led to study of a set of functions which quickly were shown to have bi- 
orthogonality properties; and their use has made possible the calculation of large quantities of 
radiation data [1-2]. 1 Reference 3 contains a brief discussion of this and related BON systems, 
together with some elementary properties easily derived, including the differential equations 
(3d order) satisfied by the functions. In this connection, it has since been proven by S. Preiser that 
the relevant differential equations are the only 3d order differential equations whose solutions are 
polynomials and have the biorthogonality property [4J. 

During the past two years, further properties of these sets of functions have been studied; 
and this manuscript summarizes the most important results. The investigations have been ex- 
tended into the approximation theory appropriate to these biorthogonal sets, and beyond this, 
some results have been obtained for related ("function fitting") types of approximations which 
have been extensively used in calculations of radiation data and other applications [5, 6]. 

The work has been largely stimulated by the intrinsic beauty of the mathematics, although the 
practical and theoretical importance of these sets of functions is considerable. And because there 
appears to be a lack of practical information about biorthogonal functions, a corresponding effort 
has been made to be reasonably complete in this presentation. 

2. General Structure 

BON "Bridges." BON systems such as those on which our attention is focused, are derivable 
from, and relatable to well-known orthogonal systems. To exhibit this clearly, we make use of 
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mathematical structures called pairings [7]. A pairing is an ordered triple (W, W ', (,)), with W 
and W vector spaces over the same field F, and (,) a functional on WX W which is linear in the 
first coordinate and antilinear in the second. The functional (,) is also frequently referred to as a 
Paring. . . . . 

Let {W, W \ (,)) and (ft, ft, [,]) be two pairings. A linear transformation, B : JF— »ft, will 
be called adjoint if there exists another transformation, B : ft ^ W such that 



\_7T,Bp] = (B7T,p), 



(1) 



for all functions ireCl and peW. A linear transformation. 0:fi— * CI, will be called self-adjoint 
if and only if 



[0£, r] = [Bt, *]* 



(2) 



for every pair of elements, o-, r, in ft. 

In what follows we assume that in addition to the two pairings mentioned, (W , W, (,)) and 
(ft, ft, [,]), we have a linear transformation,, Z?:ft— » W, its adjoint linear transformation, 
B : W— * ft, and a self-adjoint transformation, B : ft — > ft. We define 



T=BOB. 



(3) 



The situation is illustrated in figure 1. 

A I 

W 




Figure 1 



Clearly, if G is self-adjoint, so is T. For let 5, teW; then 

<7$, t) = (BeBs.t;) = [eBs, Bt] = [GBt, Bs]* = {BGBt,s)* = (ft,*s)*. 
Two systems of vectors {a,} C W, {%} C W will be called biorthogonal if and only if 

(in, uj) = 8ij 

A 

Now, assume that the elements {2t,-} C W have the property 



(4) 



and assume that {tt[\ C ft is BON to {rn} C ft. Further, assume that 

and let 

Tili = Hi. 



(5) 



(6) 



(7) 



2 The asterisk in (2) and subsequent expressions will indicate complex conjugate. 
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Then {ut} is BON to {in}. For 



(uu Uj) = (Tut, uj) = (BGBu u U j) 

= [677-;, 7Tj] = [iTu TTj] 



— °V 



(8) 



Note that in the calculation (8) it was not necessary to assume that 9 is self- adjoint; but this will 
be so for all cases of interest to us. 

We will refer to the structure of figure 1 as a "bridge." And we will refer to B and B as "bridge" 
operators. The terminology is not in use, to our knowledge; but it is convenient here. 

Operators, Spaces, and Norms. In what follows we shall make the following assumptions: 

(a) The operators B, B, and B are one-to-one and onto, and their domains are the whole spaces 
ft, W , and ft, respectively. 

(b) B and B are adjoint, and 6 is positive definite and self-adjoint, as already indicated. 

(c) {rn} C ft and {tti} C ft are biorthogonal, and 07r,- = 7r/ for all i. 

(d) Without loss of generality we may assume that the space ft is complete in the weak topol- 
ogy 3 cr(ft, H), for if it is not, we may consider its abstract completion. We also assume that the 
set of all finite linear combinations of the 7)-^ designated by (77-;), is dense in ft. 

(e) The {77,} are linearly independent. 



We now define a norm on ft as follows: 

\\a\\=[os-,&yi 2 . 

It can be shown that (9) does indeed define a norm on ft. Further, if 



(9) 



■%am, 



then 



01 






= 2>l 2 - 



We can thus map {tti} by means of an isometric isomorphism onto the set of all finite sequences 
with square norm, E%. Extending this mapping by continuity we find that ft is isomorphic to I 2 . 
Thus we obtain the following results: 



< 00 



(a)n={|)a,^ fkl 

M = 1=0 

(b) aeh=) ||S-|| 2 = |>,'| 2 - 

(c) The set {77,} is a Schauder basis for ft. 



This is the notation used, e.g., by Wilansky, ref. [8], p. 235. 
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(d) The spaces W, W, and ft are all isometrically isomorphic to / 2 and if cr=Bs, cr=Ocr 
s = Bcr= Ts, then 

W=IPII=<73^>i>*=m-I^={2N'} " 2 - (io) 

(e) The sets {ui}, {w ? }, and {77^} identified by (5), (6), and (7) are all Schauder bases for their 
respective spaces. 

While the operator T has been identified by (3) in terms of the (ft, ft, [, ]) pairing, this operator 
may be defined independently; in this case (10) relates the norm based on T to the norm of the 
other pairing of the bridge. 

This highly structured "bridge," illustrated in figure 1, is now seen to be an isometric isomor- 
phism between pairings. It enables us to move with confidence in developing and using BON 
systems. 

Kernels and Their Series Representations. In what follows, W and W will consist of functions 
defined on domains D and Z), respectively, of the complex plane. The pairing (, ) will be defined by 



(5,?) = J Zdwst*, (11) 



where Z is a suitable subset of D H D and dw will in general include a weight function. Likewise, 
ft and ft will be spaces of functions defined on domains A and A, respectively, of the complex 
plane. The pairing [, ] will be given by 



W,r ] =j [ 



cfe>OT*, (11') 



where T C A Pi A and daj will similarly include a weight function. 

The operators B and B are defined using a kernel function k, which is integrable on the domain 

(DUfl)X(AU A) as follows: 

[lK]fa)= J dw{z)k{z,t))i{z*) (12) 

and 

[Bo] (2)= J , gMt7)A-(z,t7*)o-(t7), (13) 

where Z' and f" are in D D D and A PI A, respectively. The operator 0:11 -> £1 is of the type 



[Qt] fo)=L<Mii')*(i,, V ')t(v'*), 



(14) 



where #, which is integrable on the domain (A U A), is symmetric in r) and 17', and T" is in A Pi A. 
From (12)-(14) it is clear that the operator T:W-^ W has the form 

Ts{z)=\ dw(y)t(z,y)s(y*) (15) 

where t(z, y) is a kernel function having the form 

t(z,y)=\ da)(iq) ! do>(V)£(z,^*)#(^ V)Mr, V*). (16) 
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Because the set {tti} is a basis in ft, we can write, for z fixed in D, 




k(z,r))=^km(ri). 

i 


(17) 


The functions {uj} already mentioned are given by 




Uj = Tuj = Bttj = I do) ( 7} ) A : (z, rj* ) 7Tj (17 ) ; 


(18) 


and it is clear from the series representation (17) for k(z, r)) that 




i 




so that /r (2, 77) must be given by 




/c(z,T7) = ]T tti(z)7ri(7)). 


(17') 


Likewise, from (16) and (17'), the series form for t(z,y) is given by 




t(z,y) = ^Ui(y)BQ7n = ^Ui(z)ui(y), 

i i 


(19) 



where 2, yeZ. 

3. The U* 9 U* BON System 

Reference BON System. We wish to consider functions it n {v\) given by 

/7? 2 -l \" 1 

ltf(l») = "S- " ( 20 ) 

together with another set, {77^(17)}, given by 

»-(i)-(^i)". (2Q,) 

where 17 is complex. These two sets of functions are BON for integration along a contour c parallel 
to the imaginary axis, and crossing the real axis between and l/V 2: 



^~ j c d VyZ^2 7r ™( 7 )^n( r ))= Snm- 



(21) 



Here 



[i^}" 



ays the role of a weight function. To see that (21) holds, we note that for n > m, 



the integrand has the form 

( n - /«]- 1 



air) 3 2i , 



which has no residue. Further, for (m — n)=j ^ 1, the integrand can be written in partial fraction 
form, 

201 



-2Tf- 



(r? 2 -iy 



i+l Zj Ci C\4-r,V 2s 



d + T?) ( 



'(l-T))*' 



where the sums begin with i = 2 because (a) the integrand is odd in 77, and (b) for large rj the inte- 
grand must be of order 77 ~ 3 . From this form it is obvious that the integral again vanishes; there are 
no residues at 'the poles at 77 = ±1. Lastly, a simple residue calculation establishes the result 
for m= n. 

The Spaces ft and ft. If, now, in accordance with our discussion in section 2 we define ft to be 
the set of all functions O satisfying 



#(77) = j? a n 7Tn(ri) 



(22) 



with 



2 KI 2<co > 



(23) 



and ft to be the space of all functions 4> satisfying 

DC 



(24) 



we obtain the following results: 

— - — I*. The magnitude of this variable 
will be unity on the equilateral hyperbola given by 

1 



cos 26-n = 



2hr 



(25) 



which is shown in figure 2. For values of 7] in the two regions designated I, to the right and left of 
7j 2 -l 



the hyperbola, 



V z 



< 1. Thus the series (22) will converge for all O and all 77 in regions I; and 




FIGURE 2. The hyperbola intersects the real axis at ± 1/V2. 
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$(77) can accordingly have no singularity in regions I. Any singularity must either be in region II 
or on the hyperbola. 

From (20') and (23) it is obvious that the series (24) for <£ will converge for all 17 in region II. 

The transformation 77' = V 1 — tj* 2 changes 0(77') into ^(17) and maps each singularity 



VT 



7) 



of O in II into a corresponding singularity of 4> in I. It is possible for ^(t/) to have singularities on 
the hyperbola. Clearly these are mapped into singularities of 0(77) at the same locations on the 
tj 2 -1 



1 hyperbola, since 



V 2 



= 1 implies that 



(t?'*)--1 = t? 2 = t?* 2 -1 

(t)'*) 2 7J 2 -1 7}* 2 ' 

Referring back to (19), we identify the BON-system kernel by 

#(t,, v) = 2 w-fa^-fo') = a 1 ~s (i + ~y ) - (26) 

The series converges for i) and 7)' both in region II. From a simple residue calculation we confirm 
that 

If <J>(t7) is given as in (22), then 

[e*](7)) = ^(T?)= J a„e7T„= J a„ir„(ij), 

tt=0 W = 

where the kernel of G is given in (26). Likewise, in agreement with (9) and (10) we write 

11*11 = \M = { f AjT^Wft) P = \ 1 l«»l' J }" 2 - <27) 

^ •'c ' »=0 

It is clear that the sets {77,,}. {77,,} are Schauder bases for H, Q, respectively, and that G is 
one-to-one onto. The inverse operator, G _1 is easy to write down explicitly. 

A 

The Un and U„ Functions. We shall now complete the structure outlined in section 2. 
For every nonnegative integer k we shall define an operator B^ on fl generating a sequence 
of functions which we will call {£/{j}j =0 

These operators B^ will turn out to have adjoints (in a sense to be defined presently) and we shall 
then define a system {U%}%=o BON to {£/£}£ =0 by solving the integral equations 

B,UJi = TrH n = 0,l, 

A A A 

The spaces Wu and Wk will then be constructed around the {Un} and {11*} exactly as (1 and O 
were constructed around the {tt h } and {77,,}. We shall then obtain a BON bridge for every non- 
negative integer k. 
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To carry out this program, we identify the ''bridge" kernel function as {Tr}* A >- T J** + l*l}, which 
we may also write in the form e^(— ) A ( — e _r, *M: and we apply this kernel to the 77-^17). For z^O. 

1 fco + i* 9 f <9 A ' 1 

[/*(*) =B,ir n = — *,- \eH-r — k e-w\iTn(T)) 

Z7TI Jco-ix (1— T) I d* J 



A- 7)~ 






(28) 



(l-r^V-l 

where < Co < 1/V2. The Fourier integral may be evaluated as follows: 
/a \2» + fr ] fix o 

= ( - ) " + "'2^! VI(I;) 2 " + A {|2| " + 1/2/ ^" + ,/2{|2|)} ' (28,) 

where /£» + 1/2 is the Hankel function. To this point we have not defined Un(z) at 2 = 0. An examina- 

tion of the structure of (28') reveals the reason. The operation I — I can be applied to the quantity 

in braces in (28) without difficulty, due to an interesting property of the functions z n+ ll2 K n + 1/2(2): 
they lack the first n odd powers in the MacLauren expansion. But the next derivative (k=l) not 
only generates a function odd in 2, it produces a discontinuity at z = Q. The second derivative 
(A =2) produces an even function; and the derivative of the discontinuity generates a delta function 
at 2 = 0. Higher derivatives produce functions of parity k and singular features of the same and 
higher orders at z = 0. 

The singular features at 2=0 are not, in general, characteristic of the functions to be repre- 
sented by the {/£; and they contribute to none of the relevant integrals over the £/£. Hence we 
remove them by the simple expedient of defining 

J7S0)=Bm Ufrh (29) 

Equation (41) gives these values for k = 0; for fc = l, £/ ; }(0)= 1 for all n. 

A 

We now turn to the problem of determining the functions C/J which are BON to the U%. To 

A A A 

derive these we solve the integral equation BkU$= n^, where Bk is given by 

|J"" dz{z k e-M} {« W (jfp*)*« H " W *} U«(*y=*«(v)- 
This expression simplifies to the complex conjugate form 

J ^ dur-»{*fifr)} = (^)" ±- v (30) 

From the binomial expansion of the term on the right it is clear that 

Finally, we write the BON property of these two sets of functions explicitly to exhibit the 
pairing integral, and the weight function, clearly: 
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\ j CC dz{z k e-^}U k n (z)UU2)= r dzz k e -*U k (z)(JUz) = 8 nm . (31) 

The Spaces W k and W k . The spaces Wk and Wk can now be constructed as outlined in the 
previous sections. We define 

F k (z) e W k H F k (z) = ^cti f/f; 2 H 2 < °° 
r=o 

and 

F k (z) e r fc <=) P k (z)= J <*#*; V |a*| 2 < °°. 

t=0 

It is clear that the elements F k (z) of the space W k will consist of the functions B k Q>, where Oefl. 
That is, for 2^0, 

F k (z)=(-)^(£) k JL £ A,e- ^ *(r,). (32 ) 

Note that ^(2) has parity ( — ) k . If lim Fa-(z) exists, we define F k {0) to have this value. 
The function Fk(z), which is given by 

F k (z)=% a„U*(z), (33) 

H=0 

exists and is bounded for all k, all 4> e 12 and all z/0. This is partly because of (23) and partly 
because of the fact yet to be shown that 

i£/*(z)i=o(«-n 



The inverse operator, B k l , is easy to write down and can be useful in identifying elements 

of W k : 

B^F k (z) = j*J* {^f 1 eV j e ~ lzlF ^ = *fa)- (34) 

We recall that <£>(r)) has singularities in region I, and that the power series which comprise the 
elements of H can quite generally represent functions with singularities in this region or on the 
hyperbola. Accordingly, we identify the elements of Wk by parity, by the feasibility of performing 
k integrations by parts on e~^F k {z), and by the location of the singularities of the Fourier transform 
after this modification of the integrand. 

To identify the space Wk, we first observe that we can write the operator B k in the form of a 
Laplace transformation. Referring to (30), we should have, for each F k e Wk, 



X dze-^ z z k F k (z) = - T1 7 <f>(i7), 
Jo V 



(35) 



or, eliminating the factor 17** by successive integrations by parts, 

f dze^ z (^)W*(z) = <!>(*?). (35') 
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If the function F k has even parity, and if 



iP p ^ 



has a Laplace transform whose singularities are confined to region II, then surely F k e W k . 
Further, since the operators B k , 0, andB* are all clearly one-to-one, the function F k = (B k Q B k )F k 
is then an element of W k . 

The Norm and BON -System Kernel for W k and W k . We define a norm on both W k and W k as 
already indicated: 

U^l 1/2 r x 11/2 

dz{z k e-*}F k (z)Pt(z)\ = £ \a*\*\ ■ (36) 

o J L n=0 J 

It is obvious that the {£/£} and {£/£} constitute Schauder bases for W k and Wk, respectively. 
The operator T k , with the property T k F k = F k , together with the BON-system kernel t k (z, y) 
are easily identified. If z # 0, we have 

)(i-V 2 ) 



X { J o " dyy*( -) fr (l^""'*'} **W. 

= | o °° dy(y*e-«') {eM+W (^" | Ko(Vz* + /)} £ ft (y). 
From this it is clear that for z t^ 0, 

7y^ = | [" <We" l! ") {el«W»l (J^)" | /SC (V?+7)}^(y), (37) 

where K is the Hankel function, and the term in braces is the kernel function. 

The behavior of the kernel for higher k values can to some extent be inferred from the k= 1 
case: 

(|r 2 ) I *(v7+?>- £ ^ *(V?T?), ( 37') 

which, for small values of the argument, is proportional to yzl {y 2 + z 2 ) 2 . 

4. Properties and Expansions 

Recursion Relations which Connect Systems with Different k Values. If we differentiate (28') 
with respect to z we can show that for 2#0, 

j {e-\*\U^}=-e-^Ul+\z). 

(38) 

Analogously, inspection of (30') shows that 

P <W%zO = z* +1 frS +1 (z). (39) 
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A more elaborate recursion relation, which likewise connects systems with different values 
of k, can be obtained from (28) by integration by parts. We write 4 



«*e-M0K 



IT J- 



dse 



(1 + S 2 ) W+1 J 

Integration by parts k times gives 

Carrying out one of the indicated differentiations, we obtain 



(40) 



z*e-Wf/*(z) = 



77 J- 



cbe"^ ( — I Ik 



«2n+(&-l) -2(w+l)+(fr-l) 



2n+(*-l) »2(n+l)+(fe-l) 



+ 2/7 



(i + s 2 )" +i (i+5 2 y +2 

-z^ 1 eH-l{^^i( 2 )-2^{(2) + 24^;- , (z)-f/^{U)l} (40') 

Some recursion relations which connect different functions for the same value of k will be 
included in a report to follow. 

The Trend of the U„ for Large Values of n. We now wish to examine lim U%z). We begin by 
considering the case rc = 0, for which we write 



|e-l*lf/8(2 



-r 

TT Jo 



ds 



s 2 \» 1 

1+S 2 / 1+S 2 



f x d5 (-*-Y -J-^-1- Bs±il B c/o (0) 



(41) 



Next, we change the variable of integration in (40) to // = s/vn + 1 and rewrite in the following 



form: 



z*e-'*lE/ff(z)= ; ] | | X du cos (hz V/? + 1 ) f— ) k u k - 2 {e- 1 /" 2 + R (n, «)} (42) 



where /?(rc, //)= exp {— {n+ 1) log[l + a~ 2 (/i + 1) _1 ]} — exp (— l/u 2 ). The R(n, u) gives a correction 
with trend approaching proportionality to (n+ 1) _3/2 for large n. 

We now observe that the main part of the integrand has the form 



/ d\ k 9 9 * k 

\T) {uk ~ 2e ~ 1 ' u ) =e ~ l! u ^a ki u- 2 - 2 K 



where Ik is finite and a/co is zero for k > 1. Then 



2 f °° / / rf \ * 9 

da cos (us Vn+1)( — J {u k - 2 e- l ' u ) 



du 



— I du\(. 

77 Jo 



^ aktu- 2 - 2 ^- 1 ^ 2 



:os (az V rc+ 1 )| 
where Cfc is a positive constant which, for k = 0, has the value 2 V77, in agreement with (41). 

4 This integrand has a nonconvergent component which contributes only to the (irrelevant) singularities at z = 0, and which does not invalidate the results of this 
argument. The integral on the right of (40) lacks this component. 
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We conclude that 

IW«|<^=W-*«W. ,43) 

and it is evident from this that for large values of n, for z¥=0, and for any A;, 

0S(z)=O(n-»/»). (43') 

Generating Functions. Useful generating functions for both U% and U% are easily derived- 
First, consider the sum 

,|»"k'(T^)"r^!=wf 2 W- <«> 

The singularities for this function are located at 7]= 7= =, and will be located in regions I for 

VI — u 

\u\ < 1. Fourier inversion of both sides gives (for z ^ 0) 

/ 3 \fc I HL_ 

-f - ? 4=e VT^=V tt «e-W£«W. (45) 

This may be written, 

^ -== e VI ^ = £ u»c-W£^(z); (45') 

\|z|/ \Vl- tt / ~ 

and it is clear that for \u\ < 1, this function is an element of W^. 

The convergence of (45') for |a|= 1 is important for later developments; and we establish 
this next. Note that the function in (43') is not in W k if \u\ = 1. For z # we write 

n2 \ n 



N I fioo 271 k N ( T) 2 \ 



The first term in braces gives the result of (45'). The nature of the second term becomes clear if 
we write the factor outside the braces as follows: 



2r) k _ 1 



l-7) 2 (l-u) \-u 






The term in brackets reduces to a finite sum of terms 7] k ~ 2i l{\ — u)\ plus a proper fraction. The 
Fourier transformation thus yields a finite sum V e~'*l(l — M)"*f/J+f*(^)» P ms a convolution over 

£/ v+1 (z) or ^iv+i( 2 )- All of these terms tend to zero as N becomes large, because of (43). Hence (45') 
holds, for z 7^ and for \u\ = 1, so long as m # 1. 
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To obtain a generating function for the £/£, we consider the sum 

t 2 (l - u) + u It? + 1 a + r? 2 (l - a) J 



(46) 



It is obvious that for \u\ < 1, and 1*1 < —7= (f real), the singularities of this function all lie in 

V2 

region II. Laplace inversion gives the result 5 
1 



\ t * e -t\z\ + cos 2 -J- ^ J- sin 2 a/- 

\-u yl~u yl — u Vl — u 



t 2 (l-u) + u 

= £(-')* i,u»z*U*(z). (47) 

A" = n=0 

As special cases we obtain (see footnote 5) 

-±- cos z Jr^-— f ^8(2), (48) 

1-u Vl-zz ^ 



-J- sin 2 J T JL -= J xl T — u»zW n (z); (48') 

1-w VI — h ,jto V 1 — a 

and these functions, for \u\ < 1, are elements of # () and JFi, respectively. In combination, (48) 
and (48') give 

T^ e ~ J u»fc(z)+i J^-zUi(z)\ (48") 

TAe U„ Polynomials and the Laguerre Polynomials. If in eq (45') we write u=(p— l)/p, 
the convergence of the series will be maintained for Re p ^ (1/2). Term-by-term Laplace inversion 
of this expression multiplied by p~ l is possible, and yields the result 

where L n and Hek are the nth Laguerre and kih Hermite polynomials, respectively. 

Note that (49) is the kernel for a transformation Z?a which links the spaces Wk- and Wk with 
the Hilbert space spanned by the Laguerre polynomials. 

The UJi Polynomials and Tchebycheff Polynomials. In precisely the same fashion we can 
establish bridges with other Hilbert spaces. In particular, let us write u = e ie in (45'), which is 
permissible in view of the convergence argument of section 12. The real part of the result can be 
written in the following form, which defines a function ka (z, 8): 



5 These results, in a different series form, have been given by Preiser (ref. [4]). 
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K k (z, 6) = 2 cos »«c" w WW 

- / £+lReexp i z ' e ^-^ + i(A;+l)(-7r-fl)/4[- (50) 

(V2 sin 02) I V2 sin 0/2 J 



A variant of this result is perhaps more useful in applications. Recalling that cos (nO) is the 
Tchebycheff polynomial T(cos 0), we write x— ( z J = cos 2 (612) to obtain 

r w (2*-l) = eos n[2 cos" 1 Vx]. (51) 

This polynomial system is actually one of the Jacobi family, and ve n T n (2x — 1) is orthonormal on 

the interval < x < 1 with respect to the weight function ; . Writing 

irV*(l-*) 

u = exp i[2 cos -1 V%] (52) 

in (45'), we obtain 



Kk 



(z, 2 cos- 1 Vx) = J r w (2x-l)e-l^lf/ A n U). (53) 



This transformation kernel links Wh and Wk- with the space spanned by these Jacobi polynomials. 

The UJfs and £/ie Legendre Polynomials. A bridge for the Legendre polynomials is likewise 
useful and easy to write down. A generating function for the Legendre polynomials is 

— — 1 - = V s»P n (cos 6). (54) 

{ 1 — 2s cos 6 + s-y' 2 ' ^ 

The expression on the left has singularities for \s\ = 1, so that the series converges uniformly for 
\s\ < a < 1, where a is constant. 

We write yfl{rf — 1) in place of s in (54) to obtain 



For cos 7^1, this function is in (1 (because V Pn( cos 0) converges) and we calculate its 
Fourier inverse to obtain w=0 

»=» 7rJo 1+a /i 77 a M ZZ7 « 2 t 

Unlike the preceding kernels, this integral does not appear to yield an elementary quadrature in 
terms of known functions. 

Existence of Functions with Zero Alternate Moments. One question which is implicitly an- 
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swered by the preceding analysis has to do with the existence of functions F(z), defined on 
^ z ^ oo, which are not identically zero but whose alternate moments F n are all zero beyond 
some value n — k — 2. 

Let an even function F(z) have a Fourier transform <£>{r)) with singularities only at ^ = ±7/0, 
Re 170 7^ 0. Then one can choose a number a such that F(az) has transform singularities in region I, 
and hence can be represented by a U% series. If all even moments on ^ z ^ cc G f F(z) should be 
zero, all coefficients of this series must be zero also, so that F(az) must be identically zero. Thus 
if F(z) is not zero identically, these even moments cannot all vanish. 

This argument can be extended to odd functions and their odd moments, and to functions of 
either parity which might have alternate moments zero beyond n = (k — 2), by use of expansions 
with higher values of k. The argument then applies to some component of any function whose 
Fourier transform has a singularity off the imaginary axis, and we conclude that only functions 
with singularities restricted to the imaginary axis, or at °°, can have zero moments beyond some 
finite value n for the moment index. 

Gamma Ray and Neutron Penetration Distributions. A series of investigations of the deep 
penetration trends of these distributions has shown that the Fourier transform can have singu- 
larities of the following types 

1 _c_ Vc 

(fJLo-p)*' € "-'' e ~V^' ( 56 ) 

in the important cases of isotropic or plane monodirectional sources. Here, /x , /n, /x w , c, and £ 
are all constants. By measuring distances in units of /xd 1 , /x -1 , or /x-m\ all three singularities may 
be located at p=l; and in each case the only other possible singularities are located on the real 
axis, with p > 1. 

For cases of plane, slant sources, the location and nature of the singularities is not so com- 
pletely specifiable. If distances are measured in fiy,?, where fim is the smallest attenuation co- 
efficient, we can say that the singularities are again on the real axis, with p^ 1; but we cannot 
specify them in greater detail. 

In all these cases, the even or odd component distribution functions belong to Wu according 
to the value of (— ) A , so that £/£(z) series representations should have square summable coefficients, 
and should converge to the distribution. Most approximations of this type have implicitly used 
the norms of the Un systems; since convergence of the polynomial coefficients has been used 
as a criterion of validity of the approximation. 
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